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The Generalized Langevin Equation for the Cooperative Dynamics of interacting polymer chains 
' (M. Guenza, J .Chem.Phys. v. 110, 7574 (1999)) is implemented to investigate the anomalous dynam- 

»' I ics of unentangled polymer melts. The proposed equation of motion formally relates the anomalous 

center-of-mass diffusion, as observed in computer simulations and experiments, to the nature of 
^ ' the effective intermolecular mean-force potential. An analytical Gaussian-core form of the potential 

O . between the centers of mass of two polymers is derived, which agrees with computer simulations and 

■ allows the analytical solution of the equation of motion. The calculated center-of-mass dynamics is 

characterized by an initial subdiffusive regime that persists for the spatial range of the intermolecu- 
lar mean-force potential, and for time intervals shorter than the first intramolecular relaxation time, 
in agreement with experiments and computer simulations of unentangled polymer melt dynamics. 
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The dynamics of simple fluids (e.g. monoatomic or colloidal) is driven by the fluid-mediated interactions between 
the particles as they move under the influence of the "effective" potential of mean- force. The potential of mean- force 
is related to the "static" radial distribution function g{r), which describes the probability of finding another particle 
at some distance r from the first "tagged" particle. The structure of the fluid, and the shape of g{r), depends on the 
I strength and shape of the "bare" interparticle potential, so that the nature of the interparticle interactions completely 
^ • defines the static and dynamical properties of the fluid. 

I ' For molecular fluids, the liquid structure and dynamics become more complex as they are driven by the interplay 

■ O between intermolecular and intramolecular interactions. |^,^ However, for high molecular weight fluids such as polymer 
melts, the description is again conventionally simplifled by assuming that the intramolecular forces dominate the 
intermolecular forces due to chain connectivity. In this case, the dynamics of a single chain is described by a Generalized 
Langevin Equation (GLE) where the surrounding fluid acts as a mean-field medium that interacts with the "tagged" 

^ • molecule only through the effective friction coefficient and the random intermolecular forces. Specifically, unentangled 
' polymer fluids are described by the Rouse equation, an "intramolecular" GLE where the memory function is neglected. 

I The Rouse equation has proved very successful in explaining several experimental findings of unentangled polymer 
• fiuid dynamics. Q It correctly predicts the scaling with degree of polymerization of the diffusion coefficient and bulk 
viscosity. Nevertheless, it appears to fail in describing the center-of-mass (c.o.m.) anomalous dynamics experimentally 
observed for t < Tuouse, with TRouse = Rl/Dnouse the first Rouse relaxation time. Here Rg — is the molecular 

' radius-of-gyrations, with N the degree of polymerization, and I the statistical segment length. The Rouse diffusion 
coefficient Duouse = ^bT / {NC)^ where ks is Boltzmann's constant, T is the absolute temperature, and C, is the 
' Rouse monomer friction coefficient. In the short-time regime {t < TRouse) the Rouse equation predicts linear- in-time 
^ purely diffusive c.o.m. dynamics, while the observed c.o.m. mean-square displacement follows a subdiffusive regime 
' (Ai?^(i) oc f^), with v w 0.75 — 0.9, which crosses-over to the free Rouse diffusion {AR'^{t) oc t) a.t t x trouse- 
A subdiffusive regime has been observed also in melts of entangled polymer chains for time shorter than the 
Q \ entanglement time, te (for t > te the polymer motion is restricted to a curvilinear diffusion) where the conditions for 
O ' Rouse dynamics - purely diffusive behavior - would seem to hold. In polymer solutions, the subdiffusive dynamics 
^ only appears above the overlap concentration, c* [^,^, suggesting that this effect is driven by intermolecular forces. 

rS _ The failure of the Rouse approach in the short-time regime is actually not surprising. The underlying assumption of 
I this model is that the dynamics of the surrounding fluid, and hence the intermolecular time correlation functions, relax 

■ ■ ' on a time scale much shorter than any of the intramolecular processes. This allows the minimization of the memory 

contribution resulting from the projection of the fluid dynamics onto the intramolecular coordinates, defined as the 
slow variables in the system, [ p^ In polymer melts and concentrated polymer solutions, this assumption obviously 
does not hold, since there is no separation of time scales between intramolecular and intermolecular relaxations. The 
characteristic time of intermolecular relaxation, Tinter , is defined as the time required for the polymer to escape from 
the influence of the intermolecular potential of mean-force. In polymer fluids the range of the potential is given by 
the extension of the "correlation hole" in the pair-distribution function, which is of the order of the molecular 
radius-of-gyration, Rg. Thus, Tinter appears to be of the same order of magnitude of the longest intramolecular 



X3 



1 



relaxation time, Tinter ~ trousb- As a consequence in the time region t < trousbi where the intermolecular correlation 
functions have not yet decayed to zero, we expect the dynamics to be substantially affected by the intermolecular 
effects, and the Rouse equation should not hold. 

The correction due to intermolecular effects becomes more relevant in melts of high molecular weight polymer chains 
since it depends on the number of interacting polymers to be found on average within the range of the potential, 
n « pP^/N. Intermolecular effects also increase with the number of molecules undergoing cooperative dynamics, 
when the system approaches its glass transition. 

While the theoretical approach to the structure of the macromolecular fluid has been successfully implemented 
through a rigorous formalism to include at the same level of detail both intramolecular and intermolecular interactions 
p^ , the rigorous, non phenomenologycal, upgrade of the dynamics appears quite challenging. Our approach is in 
this direction. We derive a GLE from the first-principles Liouville equation by projecting the dynamics [Q on a 
set of variables given by the coordinates of the molecules undergoing slow cooperative dynamics. Polymer melts 
are dynamically heterogeneous, and it is possible to observe interconverting regions of slow and fast dynamics. 
When projecting on the group of slow interacting polymers, the projection approximation, together with the mean- 
field description of the liquid of fast molecules, becomes well justified. The derived GLE describes the dynamics 
of a polymer chain inside a dynamically heterogeneous melt. The equation explicitly includes intramolecular and 
intermolecular forces in both the frequency term and the memory function. The traditional single-polymer GLE is 
recovered when the fluid is described as a mean-field continuum, i.e. the pair distribution function g(r) = 1. In 
the opposite limit of two interacting structureless particles, the equation recovers the well known GLE, extensively 
studied in the past. Thus, the equation holds in both the well-known opposite limits of dominant intramolecular 
or intermolecular interactions. 

Starting from the Cooperative Dynamics Generalized Langevin Equation (CDGLE), we investigate here the anoma- 
lous c.o.m. dynamics for an unentangled polymer chain inside a dynamically heterogeneous melt. In section II we 
present an analytical Gaussian core expression for the intermolecular c.o.m. potential between two interacting poly- 
mers which allows an analytical solution of the GLE through a normal modes analysis. Anomalous scaling exponents 
and asymptotic values for the c.o.m. dynamics are presented in Section III, together with numerical model calcula- 
tions. The comparison of the theory with computer simulation data of unentangled c.o.m. dynamics follows in Section 
IV. A brief discussion concludes the paper. 



II. CENTER-OF-MASS COOPERATIVE DYNAMICS, AND INTERMOLECULAR POTENTIAL 

We start from the CDGLE for the dynamics of a single polymer chain in a structurally homogeneous, but dynami- 
cally heterogeneous fluid. An effective segment of index a, in a polymer of index i, follows 



k<j b^a-^° 

- f E E / < • + it) , 
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where ^(r'(i)) is the intramolecular distribution function, — /3~^Zn(;[r^ (i), r'^(i)] is the intermolecular potential of 

mean-force between two sites (monomers) belonging to two different polymers, and /? = (/cbT)~^. Fa^^\t) is the 
projected random force (Q = 1 — P, where P is the projection operator p|,p^) that represents the random collision with 
the surrounding fluid. This Generalized Langevin equation is non linear, and couples the dynamics of the interacting 
chains through the intermolecular contributions included in the frequency term and in the memory function. The 
equation cannot be solved analytically, unless a set of approximations is introduced. 

For unentangled polymer melts the equation can be simplified, since the corrections due to the memory function 
contributions can be in first approximation neglected for several reasons. Because our approach derives the Langevin 
equation by projecting the Liouvillian dynamics over the extended set of basis functions (not only mframolecular 
polymer coordinates as in the derivation of the Rouse equation, but also miermolecular coordinates) the memory 
function corrections are further minimized with respect to the conventional intramolecular single chain description. |^ 
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The purely intramolecular equation reduces to the Rouse equation when memory function contributions are neglected. 

The Rouse equation is already a reasonable approximation to the single polymer dynamics, and is correct in the 
long-time regime. It is commonly accepted that when the memory function contribution is included, the intramolecular 
diffusive equation describes the cross-over to entangled polymer dynamics. Jl7| , p^ 

When the memory function contributions are neglected and we focus on the single chain center-of-mass dynamics 
for an unentangled polymer chain of index i, Eq.( |^) reduces to 

C^^%^ - ^T#7T n W - ril.m + Fil. W , (2) 

since the c.o.m. diffusion is decoupled from the internal relaxation dynamics. — P~^lng[r^ (t) , (t)] is the intermolecu- 
lar potential of mean-force between the c.o.m. belonging to two different polymers, C is an effective friction coefficient, 
and < Fc}n.{t) ' Fc^m.(O) >= 6(P^^5{t) the fluctuation-dissipation condition. 



A. Effective Gaussian-core potential between two interacting polymers 

To pursue on analytical solution of Eq. ( |^) we derive an analytical, Gaussian-core expression for the intermolecular 
potential between the c.o.m. of a pair of polymers. The modeling of polymer coils as interpenetrable soft spheres with 
a realistic Gaussian-core potential has been a subject of long-standing interest. In addition to the obvious scientific 
relevance |19f, such an approach could significantly reduce the computational time required in the simulations of 



polymer fluids. |20| For polymer solutions, a mean-field derivation of the potential for interpenetrating polymers 
was obtained by Flory and Krigbaum. pl| Their mean-field approach predicts a potential at contact that increases 
with N, in disagreement with computer simulations where the contact potential decreases with increasing TV. An 
alternative derivation through scaling theory ||2^ and renormalization group calculations recovers the correct 
scaling exponent in the asymptotic N oo regime. In a recent series of papers Hansen and coworkers presented an 



empirical form of the potential in the finite-size regime |24 , and showed, by comparison with computer simulations. 



that the Gaussian-core model gives a very good description of the interpolymer potential in a wide range of densities 
and temperatures. The same study also pointed out that interactions of higher order than the effective pair interactions 
do not give relevant contributions to the calculation of the thermodynamic properties of the fluid in the semidilute 
regime. 

We derive here a Gaussian-core pair potential for the interaction between the c.o.m. of two interacting polymers 
in a melt. The potential is calculated as a function of the c.o.m. interpolymer distance, R{t) = |R(i)|. If ra"* is the 
position of monomer a in chain i with respect to its c.o.m. coordinates, and r^'''' is the position of monomer b in chain 
j in the same coordinate system, the distance between a and b is given by r = r^"'-' -I- R(t) — ri* ' . At any instant the 
segments are statistically positioned around the c.o.m. following a Gaussian distribution ||2^ 

/ o \3/2 -3r2 

The distance between the polymers' c.o.m. is exactly obtained from the integral R = J dva^'^ira^) J (ir["'''5'(r^"''')r. 

In analogy with the previous calculation, the c.o.m. pair distribution function, g{R), is calculated by averaging over 
the monomer positions, the four-point distribution function which correlates two monomers on different polymers, 
and their c.o.m.. The four point distribution function is approximated by the product of three two-point distribution 
functions according to g{R) « / dra^^{ra^) J dr\f^'^{rl''^)g{r). For the monomer pair distribution function, we adopt 
the polymer thread-model form, which has been shown to reproduce many fundamental features of the polymer melt 
structure |l3l. 



9{r) = 1 + 



(4) 



Two characteristic lengths enter the equation: the local density screening length = 31 /{np*) with p* = p/^ the 
reduced fluid density, and the correlation hole length scale = Rg/V^. In the "thread" model the chain is composed 
of beads of vanishing thickness, while the "bare" interpolymer interactions are simple hard core repulsions accounted 
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for by the condition that the site-site pair correlation function vanishes at contact [g{d) = 0], while the intermolecular 
direct correlation function c(r) = co5{r). 

We solve the integral introducing the Fourier transform 

1 47r /■ „ e* ''(*) 

and integrating over the coordinates of the two monomers. The c.o.m. pair distribution function 
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The exponential terms are dominat for < N ^, which implies fc^^^ < ^ and A:^^^ < 1. This approximation 
simplifies the previous equation to 
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By keeping the first order term in the expansion of the logarithm and introducing the definitions of and ^c, the 
approximate solution for the interpolymer potential between two polymers of large, finite N becomes 

w{R) = ^logg{R) « - 1087r-^(p*)-^(iV)-i]e-3«V(^«^.) . (8) 

47rV7r ^/Np* 

w{R) is Gaussian and finite at all interpolymer distances, with a range of the order of the radius of gyration. At 
full polymer-polymer overlap, w(0) decreases with increasing reduced fluid density, p*, reflecting the transition from 
compressible to incompressible systems. It also decreases with increasing polymer molecular weight due to the 
increasing interpenetrability of the polymer chains. Eq.( is in qualitative agreement with the trend observed in 
computer simulations of semidilute polymer solutions. The last section of the paper shows that Eq.( ||) is in good 
agreement with data from computer trajectories provided by Grest for polymer melts at different temperature, 
density, and molecular weight. 

The derived interpolymer potential belongs to the class of bounded interaction potentials that do not diverge at the 
origin. |l9| Such potentials arise naturally as effective interactions between the c.o.m. of soft, flexible macromolecules 
since the c.o.m. of two molecules can coincide without violation of the excluded volume conditions. At complete 
polymer overlap, the N oo scaling limit of w{0) is not finite in disagreement with the observed behavior in 
semidilute solutions. However a comparison with computer simulations for the melt is not possible in this study since 
our data do not reach the — > oo scaling regime. 



B. Many-chain center-of-mass dynamics 

To include the pair intermolecular potential just derived in Eq.( |l|), several steps need to be undertaken. Since it 
is the distinct part of the Van Hove function that appears in Eq.( ??, we assume a time-dependence of the potential, 
which enters through the evolution of the square intermolecular distance between the c.o.m. of the two polymers, 
R^{t). The intermolecular term in Eq.( |l|) results from the decomposition of the n N body distribution function 
in the product of intra- and inter- molecular pair distribution functions, [ p^ In principle this high order distribution 
function is not equivalent to the product of two-body terms. However a decomposition in pair potentials is the natural 
approximation to adopt, and in our case it is supported by the quantitative agreement with the simulation data. The 
product of the pair distribution functions gives n equivalent terms in Eq.(|^). The intermolecular force for a group of 
71 w p*Vn interacting polymer chains is 

G{t)R{t) « ^^R-^[l - 1087r-2(p*)-2(iV)-i]e-3<«'(*)>/(4fl?)i?(t) . (9) 
87r-^7r ■' 

The intermolecular distance is approximated by its statistical average over the polymer configurational space, < 
R^{t) >, and in the many-chain description becomes < R?{t) >« nAR^{t) — 6D^t, where AR^{t) is the single 
chain mean-square displacement. D^^ = D Rouse /{p*^fN) is the collective diffusion coefficient derived below. The 
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approximate effective potential that corresponds to the force in Eq.( H), is found to be in good agreement with the 
time-dependent potential from computer simulations, as discussed in the last section of the paper. 

Eq.( with Eq.( reduces to a set of coupled equations of motion for the c.o.m. polymer dynamics 

C^^ = G{tML{t)-R^^{t)]+FQ('Ht) , (10) 

where R"''{t) = rT^ Y.%i ^clit) is the coordinate of the c.o.m. position for the dynamical heterogeneity. Eq.( [lO| ) 
has the advantage of including both the intramolecular and the intcrmolccular forces while keeping the mathematical 
simplicity of the Rouse equation, and it can be solved analytically by transformation into normal modes of motion. 

Because of the symmetry of the system, composed by equivalent polymer chains, the set of coupled equations of 
motion reduces to two independent equations in the relative, r^(i) — 2^^/'^[rimit) — viil(t)], and collective coordinates, 

C^-G{t)rdt) + FHt), (11) 
C^^FHt). (12) 

The relative and collective projected force correlation functions in the Markov limit obey the fluctuation-dissipation 
conditions 

< Fi{t) ■ F« (i') GCr^Kt - t')Sc.j , < Fl{t)F^{t') GCr^Sit - t')5o,-, . (13) 

Eqs.( |l^, ^2|) are exactly solved through standard procedures |^ , ^ , and the single chain c.o.m. coordinate, in the 
ensemble of molecules undergoing correlated dynamics, is obtained from the combination of the collective and relative 
contributions. 



III. CENTER-OF-MASS DYNAMICS: ASYMPTOTIC BEHAVIOR AND MODEL CALCULATIONS 

The single-chain mean-square displacement of the c.o.m. becomes 

Ai?2(i) = < (R,.„,(t) - R,.„,.(0))' >= < {v^t) - vdQ)f > +^ < (rx(0 - r^^)? > , (14) 

« < (r^(t) - rdO)f > +6DRouset/{pVN) . 

While the collective contribution is linear in time with a constant diffusion coefficient D^, the relative contribution 
exhibits a more complex behavior that depends on the strength of the intermolecular interaction, G(t). From Eq.( |^) 
G{0) is defined as the force at complete polymer overlap. 

< (r^t) - re(0))2 > = < r|(0) > (e'^^*' - l)' + GDRousee^'"^'^ e^^f*')^^' , 

where 

C(()) /■* _3<«^('')> rt 
a{t)^—^ / e dt' = C"' / G{t')dt' . (15) 

C Jo ^0 

Interpolymer interactions vanish at long time, t t^, where we define as the characteristic time in which G{t) — )■ 0. 
In general describes the time that a molecule has to travel to escape from the range of the intermolecular potential. 
In uncharged polymer fluids at temperature well above their glass transition, the range of the potential is comparable 
to and w Tjiouse- For t 3> t^, the single-chain c.o.m. mean-square displacement recovers Rouse diffusion. For 
t <^ Tfi, G{t) is a function of intermolecular distance < R^{t) > and Eq.( ^4|) has to be solved numerically through a 
self-consistent procedure. 

In a few well-defined conditions Eq.( |l^) is analytically soluble. At short time the single-chain c.o.m. mean-square 
displacement follows AR'^{t) w 6t/{f3NCo) + 0(0)^ < R'^(0) > tV(Co)- ^ G{0) ^ it recovers the conventional 
Langevin dynamics, while, if 10(0)1 > 0, it recovers a short-time superdiffusive behavior characteristic of Langevin 
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equations with a time dependent drift coefRcient. [ p6| , p7[ Since Langevin equations only describe overdamped chain 
dynamics, they do not hold in the short-time region where the dynamics is precoUisional. p^ ] In the present study, 
we focus on the subdiffusive and free-diffusive c.o.m. dynamics in the postcollisional regime. 

The onset of the anomalous subdiffusive regime takes place at the characteristic time ti = 2G(o) ■ This suggests that 
a system with intermolecular interaction of increasing strength "freezes" its dynamics at shorter and shorter times. 
Equivalently we can say that the motion begins to be suppressed on the largest length scale (smallest index modes of 
motion), and by increasing G(0) the dynamics is progressively frozen at shorter and shorter length scales, or higher 
and higher index modes. ||2^ The increase of the monomer relative friction coefficient has an opposite effect, since it 
slows down the dynamics and enforces the overdamped regime. 

For t 3> Ti the dynamics are different depending on the nature of the long-range potential. Since the c.o.m. in- 
terpolymer potential results from the combination of monomer-monomer interactions, and in polymer melts excluded 
volume interactions are dominant, this potential is repulsive. An effective intermolecular long-range attractive compo- 
nent can appear, however, when the single-polymer is constrained, for example because of the presence of chain-chain 
uncrossability or entanglements. In those cases, the molecule has to overcome an effective potential barrier in order 
to diffuse. Attractive effective potential are also presents in polymer blends, self-assembling polymer systems, and in 
droplets of polymer liquids. 

If we approximate the infinitesimal mean-square-displacement as linear in time, the interpolymer distance < 
Ai?2(i) >~ QDt. When introduced in Eq.( H) it yields 

2G(0)i?2 / -ffi 



a{t) = [1 e j . (16) 

For short, highly mobile polymer melts, 9Dt 3> 2i?^, the dynamics is simply diffusive at any time-scale. 

In a melt of slowly moving, long polymer chains, 9Dt ^ 2i?g. If the potential is repulsive, the single-chain dynamics 
crosses over from the subdiffusive regime, which describes the slowing down of the single-chain dynamics due to the 
presence of the surrounding molecules ("cage effect"), to free single-chain diffusion at t ^ 2R^/ {9D), where a{t) 
becomes constant, a{t) w —2G{0)Rg/{9D(). A decrease in the polymer molecular weight, or an increase in the chain 
mobility decreases the duration of the subdiffusive regime. 

The c.o.m. dynamics of slowly moving, long polymer chains interacting through an attractive long-range potential, 
follows a subdiffusive regime at intermediate time. In strongly interacting systems the relative mean-square displace- 
ment is constant, < {r^{t) — r^(0))^ >=< r|(0) > +3/[N(3G{0)], and the single-chain dynamics is driven by the 
collective contribution. The surrounding molecules act as a cage of the tagged molecule dynamics. At longer times, 
when t ~ 2Rg/ {9D), the single-chain exhibits free Fickian dynamics. This regime corresponds to the molecule's excape 
from the surrounding "cage", as it overcomes an effective energetic potential barrier. If the time scale of observation 
is shorter than this dissociation time, the c.o.m. interdiffusion is frozen, and the system behaves non-ergodically. 

The asymptotic behaviors just discussed are illustrated by numerical calculations in Figures 1 and 2. We investigate 
the influence of the intermolecular force, G'(t), in the scenarios of attractive and repulsive long-range interactions for 
long, slowly-moving polymer fluids, 9Dt S> 2i?g. In short highly mobile polymer chains the anomalous dynamics is 
less relevant. 

Eq.( [l0| ), when the definition of G{t) in Eq.( ^ is introduced, becomes non linear as the intermolecular potential 
depends on the interpolymer distance, R^{t), that changes while the system evolves in time. Thus the short-time 
dynamics depends on the choice of the initial interpolymer distance. Different values of the distance, however, mainly 
affect the dynamics in the precoUisional regime that is outside the range of validity of our equation. An estimate of 
the initial interpolymer distance, which is a function of the bulk properties of the fluid, can be obtained from the 
distribution in space of the polymer chains inside the volume spanned by the potential. The average number of chains 
in a sphere of radius r, as a function of time, is given by the time-dependent pair distribution 

n[r{t)]=N-^ f dv'pg[T'{t)] , (17) 



^0 

which gives for the initial mean-square intermolecular distance inside the range of the potential 

<R'^{0)>(x n[r{0)ydr K Rl . (18) 

Jo 

Since the intramolecular monomer distribution is Gaussian the chains inside the potential range are not uniformly 
distributed, and the resulting average intermolecular distance is of the order of the polymer radius-of-gyration. 
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To perform quantitative model calculations we specify the numerical value of a few more parameters. The model 
polymer chain comprises A'^ = 100 statistical segments of unitary length /, with the chain described by a freely jointed 
model consistently with Rouse theory. The effective temperature is set to fc^T = 1 in units of G(0), and we assume 
the density of the melt state, p = 1. The number of chains initially correlated inside the range of the potential is given 
by Eq.( |l^) as n « A:/?>tt\/nP p. n is also the number of chains undergoing correlated motion, defined from Eq.( |^) as 
the number of chains interacting through the mean-force potential. The effective time-independent monomer friction 
coefficient is set to C = 1. 

In Figures 1 and 2, we study the decay in time of the interaction strength, G(t), when its value at complete polymer 
overlap, G(0), is varied for a system of slowly moving, high molecular weight polymers. We use as an effective 
parameter the interaction strength, G'(0), that is a multiple of G(0) calculated from Eq.( ||). The features emerging 
from those plots are related to the onset of anomalous dynamics in the c.o.m. mean-square displacement as a function 
of time. In general an increase of the strength of the potential induces the onset of anomalous dynamics. 

G{t) is normalized by its own value calculated when < R'^{t) >= R^, i.e. G{t)\ii^ . For small |G(0)|, the decay of G{t) 
follows a single step mechanism that corresponds to a mean-square displacement approaching the linear single-chain 
diffusion. For large |G(0)| the behavior depends on the sign of the potential. 

For a repulsive intermolecular potential (Figure 1) an increase of G(0), for t < Tuouse, induces the onset of a stretched 
exponential decay which crosses over to Rouse behavior at t « Tuouse- In the c.o.m. mean-square displacement, the 
presence of the interactions modifies the Rouse free diffusive behavior by inducing a subdiffusive behavior at short 
time. At i ^ rRg/{9D) the system recovers free Fickian diffusion. No many-chain correlated dynamics appears at 
long time. 

To investigate the effect of an attractive effective potential we adopt G'(0) equal to a multiple of G(0), defined by 
Eq. ( ^ , assuming a minus sign. At intermediate time the system is frozen in a metastable equilibrium configuration 
with constant strength of the interactions, and constant interpolymer distance (see Figure 2). This regime corresponds 
to an anomalous subdiffusive regime {AR^{t) oc with u 1). If the thermal energy of the system is high enough, 
/3G'(0) 0, the molecule can escape from the "caged" dynamics reestablishing the effective ergodicity of the system. 
Otherwise, the system undergoes many-chain correlated dynamics, so it behaves non-ergodically. Thus, by increasing 
the strength of the intermolecular potential the theory predicts an apparent ergodic to non-ergodic transition. 

In the Figures are also illustrated the two limiting cases of a single-chain free dynamics (G(0) — 0), and for 
the "caged" collective dynamics (G(0) oo). The comparison with those curves emphasizes how the single-chain 
dynamics in strongly interacting attractive systems crosses over from the short-time single-chain dynamics to the 
collective many-chain diffusion at long time. 

IV. COMPARISON WITH SIMULATIONS 

To test our approach we compare its predictions with computer simulations of unentangled polyethylene (PE) 
dynamics in the melt state and at decreasing temperature, while approaching the glass transition (th e entanglement 
degree of polymerization = 136). PE has been extensively investigated experimentally [^9|j30[ because of its 
industrial applications, and has been chosen as a model system for computer simulations due to the simplicity 
of its monomer structure. We first compare with data from computer simulations performed by Grest and coworkers. 
[01 The simulations calculate the dynamics at constant volume and constant temperature, using the experimental 
densities at atmospheric pressure, as reported in Table I. 

We use as an input to the theory the density, polymer molecular weight, number and length of the statistical seg- 
ments from computer simulations (see Table I). Figure 3 shows very good agreement between theory and simulations, 
when one fitting parameter, i.e. the potential at complete interpolymer overlap z«'(0), is optimized. Table I reports 
the ratio between the optimized analytical values of w'{0) and the theoretical values, u'(O), calculated from Eq.( H). 
w'(0) appear to be fairly close to w{0). The agreement with the computer simulations is good for all the samples 
considered, which span a wide range of temperatures and molecular weights. 

A second approximation, implicitly introduced in the solution of the CDGLE, is the averaging of the instantaneous 
interpolymer distance in the effective potential, and the self-consistent solution of Eq.( |l^). In Figures 4a, 4b, and 
4c we present a comparison between Eq.( ^ and the time-dependent intermolecular potential between the c.o.m. 
of two polymers, extracted from computer simulations as —\ng[R(t)] for the G30HQ2 sample. The approximate 
static interpolymer potential appears to reproduce quite well the data at any time scale. This is because the time 
dependence of the potential is quite weak initially, so that the data during the first 30 ps almost superimpose, and are 
well reproduced by the same analytical formula. In Figure 4d we compare Eq.( ^ with the averaged time-dependent 
Gaussian potential, corresponding to Eq.( calculated self-consistently with Eq.( fol). The two potentials show 
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very good agreement. For time intervals larger than 30 ps, when the potential starts to differ significatively at short- 
distance, the molecules have already diffused to large intermolecular distance and the approximated analytical form 
of the potential still reproduces well the data. 

From the analytical solution of Eq.( |l^ we calculate the c.o.m. mean-square displacement. The overlap value 
of the intermolecular potential, w'(0), is fixed from the previous analysis. From the long time diffusion coefficient 
we derive the value of the effective friction coefficient, reported in Table I. Figure 5 shows that the comparison 
between simulations and analytical theory gives an excellent agreement in the whole range of space and time scales 
investigated. This suggests that both the self-consistent procedure and the approximation of neglecting the memory 
functions contribution are reasonable first-order approximations for unentangled melt dynamics. 

A. Unentangled undercooled polymer melt center-of-mass dynamics 

We now compare our approach with data from computer simulations of undercooled polymer systems. In a recent 
paper Binder and coworkers [|T| investigated by computer simulations the onset of short-time anomalous dynamics for 
a melt of unentangled polymer chains. They performed Dynamical Monte Carlo computer simulations at decreasing 
system temperature, using the bond fluctuation model of a chain of = 10 statistical segments. For this model, 
which is found to reproduce correctly the dynamics of polymer melts. Binder reports the unit of the single time step (1 

Monte Carlo Step « 10^^'^ s) and the space unit (1 lattice unit w 2.3 A)- The chain of = 10 effective segments 

o . 

simulated here, corresponds to a polyethylene chain of iV « 50 bonds of length I — 1.54 A- |p4| Since the model used 
to simulate the "freezing" of the dynamics privileges stretched bond configuration, the polymer radius-of-gyration 
slightly increases with decreasing temperature. 

We fit the c.o.m. mean-square displacement having as an input to the equation the simulation parameters. The 
monomer friction coefficient is obtained from the diffusion coefficient in the long-time Rouse regime and is kept 
constant in all the calculations. We do not include the memory function corrections |Q which in our calculations 
are negligible for N < N^. Because of the approximations involved in the calculations, we don't strongly rely on 
the precision of the numerical values obtained from the fitting procedure, however the theory makes some interesting 
suggestions on the physics involved in the process. 

In Figure 6 we show that the fitting procedure gives a reasonable agreement between the theory and the data in 
the complete range of time (length) scale investigated. At high temperature (/csT = 1,0.4) the dynamics is driven 
by repulsive interactions in agreement with our previous calculations. The repulsive potential induces a subdiffusive 
regime that crosses over to the single chain free diffusion. In the simulation the decrease of the diffusion coefficient, 
follows the empirical Vogel-Fulcher law ||3^, while the Rouse equation holds for the monomer local dynamics. To 
obtain a good quantitative agreement with the simulations we use two fitting parameters: the strength of the potential 
at complete polymer superposition, w{0), and the number of polymer chains undergoing cooperative dynamics, n. The 
best fit with the data is found for n = 13 and n = 15, respectively at ksT = 1 and ksT = .4. These values are very 
close to the theoretically predicted value. For temperatures close to the melt state the system behaves consistently 
with the study presented above. 

When the temperature of the system is further decreased {ksT = 0.23,0.19) the subdiffusive behavior appears 
to have a different physical origin than in the high temperature limit. The simulations show an "extremely slow 
motion", visible also in the monomer mean-square displacement, which is defined to be a "cage effect". |3^ From 
our calculations this regime emerges as a diffusive process involving the correlated dynamics of many-chains. In this 
regime the polymer chain is caged by the presence of the other polymers and cooperative many-chain dynamics takes 
place, which indicates an attractive effective potential. While the monomer friction coefficient is constant from the 
high temperature regime, the number of polymer chains undergoing correlated dynamics increases with decreasing 
system temperature, in agreement with experimental studies. |l^j3^,Q This effect would indicate an increase of 
the range of the dynamical correlation that would correspond to longer range oscillations in the "real" mean- force 
potential and in g{r). From the fit to the long-time collective diffusion coefficient we obtain n « 25 for ksT = .23 
and n w 300 for fc^T = .19. These results suggest that the freezing of the dynamics is related to a rapid increase in 
the number of polymer chains involved in the correlated dynamics. In this region of temperature the polymer has to 
escape from the influence of the intermolecular potential, overcoming an effective potential barrier, through a process 
which is highly cooperative, and this becomes increasingly difficult at lower temperature. 
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V. CONCLUSIONS 



The presence of effective intermolecular interactions significantly modifies the single-chain dynamics in the short 
time, local space regime, To ^ i ^ trousc- The onset of anomalous dynamics, observed in computer simulations of 
polymer fluids in the short-time region where the Rouse equation breaks down, appears to be quantitatively related to 
the presence of the effective mean-force potential. The spatial extent of the anomalous dynamics corresponds to the 
range of the mean- force potential (w Rg in polymer melts well above Tg), while the time scale that characterizes the 
anomalous dynamics is given by the time necessary for a polymer to escape from the influence of the intermolecular 
potential {t w TRouse)- 

The subdiffusive dynamics at time scale ri <C t ^ TRouse , corresponds to a slowing down or arrest of the dynamics 
due to the presence of intermolecular forces. It appears to be a cross-over regime to the long-time diffusive dynamics, 
consistent with computer simulations where a wide range of exponents is measured . 

The sign of the intermolecular potential qualitatively modifies the long-range anomalous dynamics. The theory 
quantitatively reproduces the experimental data in the complete range of time and spatial scales investigated. 
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TABLE CAPTION 



Table I: Simulation and fitting parameters. 



Polymer 


T[K] 


p [g/cm^] 


l[A] 


w(0)/u;'(0) 


Initialslope 


C [10^ 


dyn s/cm] 


C10H22 


298 


0.7250 


3.26 


2 


0.97 




0.25 


Cl6-ff34 


298 


0.7703 


3.84 


.4 


0.94 




0.44 




323 


0.7531 


3.75 


.6 


0.96 




0.30 




373 


0.7187 


3.66 


.8 


0.97 




0.18 


C30HG2 


400 


0.7421 


4.02 


.3 


0.90 




0.21 




400 


0.7570 


4.18 


.3 


0.85 




0.27 



FIGURE CAPTIONS. 



Figure 1: a) Normalized intermolecular constant versus normalized time, for various strengths of the repulsive 
interaction, G"(0). G'(0) is the effective strength at complete polymer overlap from Eq.( |^). From the top curve 
(close to the single-chain diffusion) to the bottom curve: G"(0) = 0.05 G{0), 0.1 G(0), 0.2 G(0), G(0). G(0) is the 
intermolecular interaction at initial conditions: fcsT = 1 in units of G(0), N — 100, n = ^n^/Nl^p, p — 3 and 
_R^(0) « Also shown is G{t) for the chain undergoing single chain diffusion, b) Normalized single chain center- 
of-mass mean-square displacement vs normalized time, at various strengths of G'(0), same values that in Figure 1 
a). ^ 

Figure 2: a) Normalized intermolecular constant versus normalized time, at decreasing strength of the attractive 
interaction, G'(0). From the top curve (close to the single-chain diffusion) to the bottom curve: G'(0) = —0.05 G(0), 
-0.067 G(0), -0.08 G(0), -0.09 G(0), -0.1 G(0), -0.2 G(0), -G(0). Same initial conditions than in Figure 1. Also 
shown is G{t) for the chain undergoing collective dynamics, and for the single chain diffusion, b) Normalized single 
chain center-of-mass mean-square displacement vs normalized time, at decreasing |G'(0)|, same values than in Figure 
2 a). 

Figure 3: Comparison between the analytical expression of the soft-core center-of-mass Gaussian potential, Eq.( ||), 
and computer simulations data, vs. the c.o.m. interpolymer distance normalized by the polymer end-to-end distance, 
i?ete, from computer simulations. 

Figure 4: Comparison between the analytical equation of the soft-core center-of-mass Gaussian potential, Eq.( ||), 
and the distinct part of the Van Hove distribution function, — liig[R{t)], from computer simulations of Cs,oHq2. Plot 
vs. the c.o.m. interpolymer distance normalized by the polymer end-to-end distance, Rete, from computer simulations. 
Curves at a) 10 ps, b) 20 ps, and c) 30 ps. d) Comparison between Eq.( ||) and the time-dependent Gaussian potential 
used in Eq. ([To|) , calculated self-consistently, for 10 ps (square), 20 ps (triangle up), 30 ps (triangle down). 

FIG. 5: Center-of-mass mean-square displacement as a function of time. Best fit of the molecular dynamics 
simulation data (filled circle) with the Rouse equation (dashed line), and with the intermolecular diffusion equation, 



Eq.( 10), for unentangled polymer melts (full hue). The short-dashed hnes indicate the longest Rouse relaxation time, 

^Rouse • 
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Figure 6: Center-of-mass mean-square displacement as a function of time at fcsT = 1 (diamond), ksT = 0.4 
(triangle), ksT = 0.23 (square), ksT = 0.19 (circle). Best fit of the simulation data (ref. |3l)) with the Rouse 
equation including intermolecular interactions. 



11 





Figure 1 Guenza 




Figure 2 Guenza 




*^16'^34 

T=298K 




^30 "^62 

T=400K 




^161^34 

T=323K 



P u 

^44*^90 

T=400K 




Figure 3 Guenza 




Figure 5a Guenza 




Figure 5b Guenza 




4^ \^ I I I I I I 

-2 2 4 

log(t) 



Figure 6 Guenza 



